arXivil509.01049v2 [math-ph] 16 Dec 2016 


The volume of Gaussian states by information geometry 


Domenico Felice/’[^Ha Quang Minh,^’|^ and Stefano Mancini^ ’[^ 

^School of Science and Technology, University of Camerino, 1-62032 Camerino, Italy 
INFN-Sezione di Perugia, Via A. Pascoli, 1-06123 Perugia, Italy 
^Istituto Italiano di Tecnologia 
Via Morego 30, 1-16163 Genova , Italy 

We formulate the problem of determining the volume of the set of Gaussian physical states 
in the framework of information geometry. This is done by considering phase space probability 
distributions parametrized by their covariances and endowing the resulting statistical manifold with 
the Fisher-Rao metric. We then evaluate the volume of classical, quantum and quantum entangled 
states for two-mode systems, showing chains of strict inclusions. 
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I. INTRODUCTION 

States of physical systems in classical and quantum mechanics are represented by very different mathematical struc¬ 
tures, nevertheless analogies appear at certain points of comparison [T]. Classical states are depicted as probability 
density functions (pdfs) in phase space, whereas quantum states are described by density operators defined on Hilbert 
spaces [2]. In fact, the notion of phase space is often gotten rid of in quantum mechanics because of the non¬ 
commutativity of canonical variables. Nevertheless, phase space can be considered as common playground for both 
classical and quantum states when one employs for the latter a description in terms of the so-called quasi-probability 
distribution functions, such as the Wigner function [ 3 . Then, one can address the computation of the volume of 
different classes of states in the phase space framework. The issue of the volume of sets of states is of uppermost 
importance. It can help in distinguishing classical from quantum states as well as to find separable states within all 
quantum states. Separable states are the states of a composite system that can be written as convex combinations 
of subsystem states, in contrast to entangled states [1]. Determining the volume of physical states is also relevant for 
defining “typical” properties of a set of states. In fact, to this end, one usually resorts to the random generation of 
states according to a suitable measure stemming from the volume of states 

Describing the geometric properties of sets of states is intimately connected with the evaluation of their volumes. 
The sets of classical and quantum states are both convex sets. In hnite dimensional systems, several metrics are 
introduced in order to compute the volume of physical states. Due to their own nature as pdfs, classical states 
can be distinguished by the well-known Fisher-Rao metric [S]. Quantum analogue can be found in the setting of 
pure states, where the Fubini-Study metric turns out to be proportional to the Fisher-Rao metric [7]. However, for 
quantum mixed states there is no single metric [8]. Several measures have been analysed, each of them arising from 
different physical motivations and advantages [9]. Such different measures have been proposed on the set of density 
matrices acting on a finite-dimensional Hilbert space; a very natural one employed the Positive Partial Transpose 
(PPT) criterion [1] to determine an upper bound for the volume of separable quantum states and hguring out that, 
for any composite quantum system, it is different from zero regardless of the number of subsystems it contains and its 
(finite) dimension [lOj . Other important measures include the Hilbert-Schmidt measure, the Bures measure and the 
measure induced by partial trace on composite systems. All of them use techniques from geometric functional analysis 
and convex geometry to estimate the volume of separable quantum states [IIHIS]. Finally, a generalization of the 
Hilbert-Schmidt measure and the volume induced by the partial trace on composite systems, the so-called a-volume, 
showed that the PPT criterion is not precise for large dimensions of the Hilbert space M- This is an evidence that 
with increasing Hilbert space dimensions, the procedure to test the separability becomes more and more difficult to 
implement. Additionally, when going to infinite dimensional systems (often referred to as “ continuous variable” -CV- 
systems), problems also arise from the non-compactness of the support of states. 
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Thus, on the one hand, we have the difhculties in analysing infinite dimensional systems, while on the other hand 
we still lack a unifying approach for evaluating volumes of classical and quantum states. To deal with these problems, 
we propose to exploit Information Geometry. This is the application of differential geometric techniques to the study 
of families of probabilities m- As such, it can be applied to Gaussian states, be they either classical or quantum. 
Indeed, Gaussian classical states are pdfs in phase space and the same is true for Gaussian quantum states, which 
are pdfs coming from Wigner functions in phase space [16]. The main reason for the focus on Gaussian states is that 
they are ubiquitous in physics, mathematics and information theory (see e.g. HD- 

Very recently, a method based on the extension of the Hilbert-Schmidt measure has been proposed HU to evaluate 
the volume of Gaussian quantum states, which is not applicable however to the classical states. In the present work, 
we exploit methods of Information Geometry in order to associate a Riemannian manifold to a generic Gaussian 
system. In such a way, we consider a volume measure as the volume of the manifold associated to a set of states of 
the system. More specifically, we start by considering N identical and indistinguishable particles, i.e. bosonic modes 
characterized by their positions and momenta and we assume that a Gaussian pdf with zero mean value describes 
the whole system state. Such a pdf is characterized by a set of parameters, i.e. the entries of the covariance matrix 
(depending on their values we can have various classes of states). Then, thanks to these parameters, to each class of 
states is associated a statistical model which turns out to be a Riemannian manifold endowed with the well-known 
Fisher-Rao metric (see also [H]). We are able to overcome the difficulty of an unbounded volume by introducing 
a regularizing function stemming from energy bounds, which acts as a form of compactification of the support of 
Gaussian states. We then proceed to consider a different regularizing function which satisfies some nice properties of 
canonical invariance. Finally, we find the volumes of classical, quantum, and quantum entangled states for two-mode 
Gaussian systems, showing chains of strict inclusions. 

The layout of the paper is as follows. In Sec. |n| we recall the phase space representation of both classical and 
quantum states. Then, in Sec. |III| we present a volume measure for Gaussian states based on information geometry. 
Sec. IV is devoted to the regularization of the introduced volume measure. Applications to bipartite states of two 
mode systems are discussed in Sec. Finally we draw our conclusions in Sec. |VI| 


II. PHASE SPACE REPRESENTATION OF STATES 

The phase space F of V identical and indistinguishable particles (i.e. bosonic modes) is the 2V-dimensional space 
of allowed real values for the canonical position and momentum variables ^ = {qi,Pi, ..., 9Ar,PAr)^ of such modes (by 
T we denote the transpose). 

A classical state for such a system of N modes is represented by a pdf in F, namely 

p:F^K+, = (1) 

As a particular case, when p becomes a Dirac delta — ^q) we have a pure state whose values of position and 

momentum variables are (deterministically) given by ^o- Throughout the paper we will consider F = and the 
integration is performed on when not otherwise specified. 

The probability density function in Q can be considered as originating from the characteristic function Xp('r), 
through the Fourier transform, 


p{0 = 


/ 


dr e *^""Xp(t), 


( 2 ) 


where i is the imaginary unit and r G 

The set of all (mixed and pure) states is a convex set, that is if Pj(^) for j = 1, 2,... represent states and {Pj}j is 
a probability vector, then 
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is still a possible (mixed) state. Only pure states cannot be decomposed in a non trivial manner as convex sum of 
other states, so they are the extremal points (or extremal elements) in the space of all states. 

The quantum analogue of pdf is the density operator p defined on the Hilbert space % = L^(IR)®^ associated to 
the V-mode system. The canonical position and momentum variables become operators qk,Pk, k = 1,..., on H 
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with the commutation relation [qkiPk\ = *• Setting R 2 k-i ■= Qk and i?2fc := Pk these relations are summarized as 


Rk,Ri 


= i flki, where flu is the kl entry of the antisymmetric 2N x 2N matrix 


Af / V 

^ = 0 ( -1 J ) ■ (3) 

This induces a symplectic structure on the phase space T, meaning that a bilinear form w : T x T —> K exists, w being 
non-degenerate and skew-symmetric. 

A phase space representation of the state p can be given by means of the Wigner function defined as in [3] 

J +yi,<7i - 2/1, ■ • ■ +2dv,<7Ar - 2/7v)rfyi ■ •-dj/AT- (4) 

Here, p{qi + yi,qi — yi, ■ ■ ■ ,qN + yN, qN — yN) is the position representation of the density operator, i.e. the represen¬ 
tation of the operator p on the eigenvectors of the operators qk, k G N}. In such a way, the Wigner function 

turns out to be defined over the 2Wdimensional phase space T and would be the analogous of classical pdfs p{^). 
Nevertheless, the Wigner function is not a pdf because it can also assume negative values. Hence, it is often called a 
quasi-probability distribution function. 

Yet, as a proper pdf, the Wigner function can be considered as originating from the characteristic function Xp(t) 
through the Fourier transform, 

W{f) = J dr Xp(r), (5) 


where 


XpiO ■■= tr 


pm 


( 6 ) 


and 


D{0 ■■= exp 


idkqk+PkPk) 

k 


A. Gaussian States 

Gaussian states are those for which the characteristic function is a Gaussian function of the phase space coordinates 
namely 


or 

where V is the 2N x 2N covariance matrix and p G the first moment vector (recall that a Gaussian state is 
completely determined by V and p). 

Although formally identical, Eqs. Q and (|^ differ by the conditions imposed on the covariance matrix. In fact, 
for classical states V is symmetric and strictly positive definite, i.e. F > 0. Yet, not all symmetric, positive definite 
matrices correspond to the covariance matrices of quantum physical states. In fact, due to the non-commutativity of 
canonical operators we have [Bun] 

Theorem II. 1. A real, symmetric 2N x 2N matrix F > 0 describes a Gaussian quantum state if and only if 


(7) 

( 8 ) 


F-hifI > 0. 


( 9 ) 
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Relation (|^ is equivalent to Schrodinger’s formulation of Heisenberg’s uncertainty principle |21j . Furthermore, 
the if and only if in Theorem m is a peculiarity of Gaussian states as can be seen e.g. by Hardy’s formulation 
of Heisenberg’s uncertainty principle [HI US]. In any case, it is trivial to show that Eq. (|^ implies the positive 
definiteness of the matrix V, whereas the converse is not true. 

The Gaussian form of the characteristic functions Q and ^ reflects on the corresponding phase space represen¬ 
tations p(^) and VF(C) by Eqs.([^ and ([^, which we can commonly write as 


P{0 


(27r)^Vd^’ 


( 10 ) 


Here we set ^ = 0 since the first moments are irrelevant for most of the physical properties of Gaussian states. Notice 
that the Wigner function, being in such a case a Gaussian function, is a true pdf. 


Among quantum states we can also distinguish between separable and entangled states @]. To this end, it would 
be helpful to employ the partial transposition. It follows from the definition of the Wigner function Q that on phase 
space, transposition corresponds to the transformation that changes the sign to all ps coordinates and leaves the qs 
unchanged 


{qi,Pi,---,<iN,PN) ^ ^iqi,Pi,---,qN,PN) ■= {qi,-pi, ■ ■ ■ ,qN,-PN) ■ (n) 

Consider now a composite Gaussian system with two subsystems A and B; let V be the covariance matrix describing 
the whole system, and Va and Vb be the ones describing subsystems A and B, respectively. Denote by := A0 id 
(resp. Ab = id 0 A) the positive partial transposition in the A’s system only (resp. B’s system only). Then a 
necessary and sufficient condition for the separability of the system is given by the following theorem [24] 

Theorem II.2. A Gaussian state described by the covariance matrix V is separable if and only if there exist covariance 
matrices Va and Vb such that 


V>Va®Vb. 


( 12 ) 


Unfortunately this theorem is not easy to verify in practice since doing so requires looking for covariance matrices 
Va and Vb satisfying (12). 

Nevertheless, if we consider a two-mode Gaussian system (N = 2), the criterion to distinguish separable from 
entangled states simplifies into [55] 

Theorem II.3. The 4 X 4 symmetric matrix V satisfying the condition ([^ describes a separable state if and only if 

V + in>0, (13) 


where V = AbVAb, with 


^B{qi,Pi,q2,P2) = {qi,Pi,q2,-P2)- 


(14) 


III. A VOLUME MEASURE BASED ON INFORMATION GEOMETRY 


Let us consider the family S of Gaussian pdf with zero mean in the 2V-dimensional phase space T. Each such pdf 
takes on the form P{£,) of Eq. (10) and may be parametrized using m < N{2N 0 1) reaLvalued variables ,... 0"* 


(the nonzero entries of the covariance matrix) so that 


s-.= lp{^) = p{^-e) = 


(27r)^y(teU^’ 


0 e 0 


(15) 


where 0 is a subset of M"* obtained by requiring some specific constraints on V{9), and the mapping 9 >->• P{.;9) is 
injective. In such a way, S turns out to be an m-dimensional statistical model (in fact a Gaussian statistical model). 
The parametrization is provided by the entries of the covariance matrix V = by defining 9^' = V^^, with 

fi-2 

I = 'y^{2N — r) + v — fi+l and 1 < I < m (the summation over r disappears when /i = 1). 

r—0 
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Given the statistical model S of Eq. (15), the mapping ip : S ^ K™ defined by ip{P{.;9)) = 9 is injective and it 
allows us to consider ip = [0^] as a coordinate system for S. In addition, we assume that a change of coordinates 
'0:0—)' ^/’(0) C M"* is such that the set \ k € where k is the set of new coordinates given by 

K* := ip{9^), represents the same family of probability functions as 5 = {P(.;0) | 9 G 0}. Moreover, we also assume 
that such a change of coordinates is differentiable. Thereby, S can be considered as a C°° differentiable manifold, 
called statistical manifold m- 

ft—2 

Remark III.l. From here on, we assume that 9^~^^ = 0 with I = — r) + 1, for p, = 2k + 1. This implies that 

r—0 

there is no correlation between position qk and momentum pk of the kth mode, for all k = 1,..., N. 

Consider now a point 0 S 0; then, the Fisher information matrix of 5 at 0 is the m x m matrix g(9) whose entries 
are given by m 


9^9)-.= f dx P{e, 9) d^inPi^-9)8,. in P{^; 9), 


(16) 
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with 8f. standing for The resulting matrix g(9) is symmetric and positive semidefinite. Yet, we assume from 

now on that g{9) is positive definite. In such a way, we can endow the parameter space 0 with a Riemannian metric, 
the Fisher-Rao metric, given by G{9) := g^i,{9)d9^ (g) d9’', with g^,v{9) as in Eq. ( [T^ . With this metric, the 

flU 

manifold Jxi •= (0, G{9)) becomes a Riemannian manifold. 

The parameter space 0 in (15) does not coincide, in general, with the whole linear space M™. The central issue is 


that by requiring that the covariance matrix V (0) satisfies some specific conditions, 0 can represent different states 
of the physical system. 


Definition III.l. Consider the Gaussian statistical model S = {Pi^',9)} in Eq. (15). Then, the classical states of 
the physical system are represented by the parameter space 0 given by 


0c,....c := {9 G M’"|E(0) > 0}. 

The quantum states, instead, are represented by means of the following parameter space 0, 


0„ 


:= {0 G K'"|E(0) > 0}. 


If the physical system is composed by two subsystems A and B, then its separable states stand for 


0„ 


:= {0 G K™|E(0) > Va© Vb}. 


Finally, in this case, the entangled states are given by 


^8 entangled * ^8 quantum 


zparable ’ 


(17) 


(18) 


(19) 


( 20 ) 


Remark III. 2. Eq. ( [l7| ) represents all the possible classical states of the physical system described by the 2N- 
dimensional phase space T. A labelling permutation a of the system’s modes acts on the pdf P{^; 0) by a permutation 
congruence of the covariance matrix: V{9) —>■ n^E(0)n, where 11 is the permutation matrix corresponding to a. Now, 
n^E(0)n is still positive definite; so, the parameter space 0cWc has a permutation invariant form. 

Eq. ( |18[ ) represents all the possible quantum states of the physical system described by the 2N-dimensional phase 
space r. It is well-known that the uncertainty relation E(0) + > 0 has a symplectic invariant form f^Sjl . i.e. given 

any symplectic matrix S, then also S'^V{9)S + > 0 holds true. Thus, the parameter space 0,„„„t„„ has also a 

symplectic invariant form in addition to the permutation one. 

In general, the definition given in Eq. (19) is not operational; indeed, to use it, it is necessary to prove the existence 


of the matrices Va and Vb • However, when dealing with two-mode systems, such a criterion becomes useful in practice 
by providing a necessary and sufficient condition to distinguish separable states among all the quantum states. 

Going on to multipartite systems, the task to describe separable states becomes harder and harder. In fact, a general 
criterion is still missing. 
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A. The Volume measure 


From Def. III.l we see that, each different set of Gaussian states is associated with a Riemannian manifold. Thus, a 
natural volume measure for a set of states is the volume of the associated manifold. 

Definition III.2. Consider a physical system of N modes in a Gaussian state. Let 0 he the parameter space as in 
(III.l) and A4 = (0,G(O)) be the Riemannian manifold associated to the class of Gaussian states 0, with G(9) being 
the Fisher-Rao metric. Then the volume of the physical states represented by 0 is 


V{V) := [ d9^ydet g{9), 

Je 


( 21 ) 


where g{9) is the real symmetric matrix with entries given by (16). 


Given the formal definition of the Fisher-Rao metric tensor (16), in order to apply it in practice, we consider a 


clearer analytical relation between the components of the metric G{9) and the covariance matrix V{9). 


Theorem III.l. The entries (16) of the Fisher-Rao metric are related to V by 


g^, = ^tr [V-i (d^V) V-^ {d^V)] 


( 22 ) 


for every G {1,..., to}. 

Such a relation is well-known in literature (see for example [Ml ET]); however, here we propose an alternative 
derivation (see Appendix [A|). 

At this point, we proceed to show some properties of the volume defined in Definition |III.2| Given Remark 
III.2[ we would require the volume in (HI) to be invariant under labelling permutations of the modes. So, consider 
a point ^ = {qi,Pi,... tPntPn) = (fi)---jCiv) € T and a permutation a : {1,...,A^} -G {1,...,A^} such that 
irr = (Ccr(i)) ■ ■ • jicT{N)) is Still a point in the phase space F. At the level of pdf in ( |Io| ), such a permutation acts by 
transforming the covariance matrix V {9) in the following way 

v'{9) = V{9) n, (23) 


where V and V are the covariance matrices of the state described by variables ^ and respectively, and 11 is the 


with Cj denoting a row vector of length 2N with 1 in the jth 


permutation matrix given by 11 = (eo.(i),. •., e£,(Ar))^, 
position and 0 everywhere else. 

Another feature we would require is the invariance of the volume measure in (21) under symplectic transformations 
S G Sp(2A^,K), i.e. S such that 


ns = n, (24) 

where n is the antisymmetric matrix defined in ([^. This requirement is motivated by the fact that from Eq. ( |24[ ) it 
follows that the uncertainty relation ([^ has an Sp(2iV, K) invariant form. 

The following proposition shows that both of these properties are satisfied by the volume V defined in Definition 

EU 

Proposition III.l. If there exists a permutation matrix 11 (resp. a symplectic matrix S) such that V' = 11^ V 11 
fresp. V' = V S), then 

V{V') = V{V). (25) 

In fact, the Fisher-Rao metric is invariant under more general transformations, namely the congruent transforma¬ 
tions defined by elements of the general linear group (see for example P5H5D] 1. However, for our purposes, it suffices 
to consider the behavior of the volume under the congruent transformations defined by permutation and symplectic 
matrices. See Appendix jB} 
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IV. REGULARIZED VOLUME 


In general, the integration space 0 given in Del. is not bounded. However, this is not the only reason for 

the possible divergence of the integral (21). Indeed, let us recall that such an integral is computed by means of the 


volume element coming from the Fisher-Rao metric G(0), which is 

VQ := ^/detg{^ d6^ A ... A 


Here, by Eq. (22), the entries of the mx m symmetric matrix g{0) can be written in the form 

g^,^ = F(y) (detU)"\ 


(26) 


(27) 


where F(y) is what is left after grouping the common factor (detU) in (22). Such a factor comes from the 
well-known relation V~^ = (detU)~ adj(U), where adj(U) denotes the adjunct of the matrix V. 

Hence, we have 


det5(6») = 


(detU(6i))" 


;FiVi9)), 


(28) 


where F{V(9)) denotes a non-rational function of the coordinates 9^, ... ,9"^' 


0 in Def 


HI.l 


Now, f rom E q. (28) it is clear that the reasons for the possible divergence of the integral in (21) are twofold: the set 
is not compact because the variables 9‘ are unbounded from above, which makes the quantity F{V{9)) 


divergent; furthermore, det 5 ( 0 ) diverges since det U approaches zero for some 9^ € 0. 

It is then necessary to introduce a regularizing function d>(V) which eliminates these possible divergences. It 
should supply a kind of compactification of the parameter space and excludes the contributions of 9‘ making det ^(0) 
divergent. Such a function might stem on physical arguments related to finiteness of energy. 

In general, for an arbitrary Gaussian state p with zero first moments, the trace of the covariance matrix is directly 

linked to the mean energy per mode, namely £ = (^) [SI]- Thereby, we propose to bound the parameter space 

with a suitable energy value of the state. To this end, we define a regularizing function as 


<P(V) ■= H(E - tr(U)) log [1 -f (det U)" 


(29) 


where ff(-) denotes the Heaviside step function and E is a positive real constant (equal to 2N£). 

With this regularizing function, we arrive at the following volume for sets of states: 

Definition IV. 1. Given a set of Gaussian states represented by a parameter space 0 as in Def. we define its 

volume, regularized by the functional d>, to be 


V^{V) := / <P{V) JZG, 


(30) 


where vq and d>{V) are given by Eqs. (26),(29) respectively. 


The integral in ( |30[ ) is now meaningful. Indeed, we have the following results. 

Theorem IV. 1. Let E denote the constant m x m matrix defined by 

E^,y = ^tr[(5^V)(5,.U)], l<p,v<m. 

Let aA]{V) denote the adjunct matrix ofV. The Eisher-Rao information matrix g satisfies 

2m / -I \ 2m 


detg< "det(£;) 


1 


det(i?). 


(31) 


(32) 


,Amin(H)^ 

where Amax[adj(lA)] denotes the largest eigenvalue o/adj(U) and Amin(U) denotes the smallest eigenvalue of V. 
Proof. See Appendix [C] □ 



















Corollary 1. The regularized volume element satisfies 


$(l/)v/d^ < H{-E - tr(y))A™,,[adj(C)] 


log[l + (det Vy 
(det F)"* 


Consequently, the integral 


$(C)\/det gdd, 


is well-defined and bounded for any measurable subset 0 C 
Proof. See Appendix [Cj 


over which V is positive definite. 


(33) 


(34) 


□ 


Remark IV.1. Setting the energy of the Gaussian states to be smaller than or equal or to E results in an upper bound 
for parameter space 0. Furthermore, the singularity occurring when det V goes to zero is eliminated by the logarithm 
rfji log(l “h 

log [1 + (det Vy^], using the well-known relation lim - = 1. So the integral (30) is now meaningful. 

a:- >0 X — 


However, the regularizing function d>{V) of Eq.(29 \ is not invariant under symplectic transformations. Indeed, 
consider S € Sp(2N,R) and V' = S'^ V S; then, tr(V^) = tr(S' S'^ V), which is in general not equal to tr(V). 


From Remark IV. 1 the following issue arises: the well-known Williamson’s Theorem [32] states that given a positive 


dehnite and symmetric 2N x 2N matrix V, there exist a symplectic 2N x 2N matrix S and a diagonal and positive 
dehned 2N x 2N matrix D such that V = S^'^DS |36|. This implies that two Gaussian states are similar under 
the congruence transformation via a symplectic matrix. The regularizing function ^{V) is not invariant under such 
transformations. 

In order to overcome this problem, we propose a different regularizing function T(V) which is devised by exploiting 
the relation given by Theorem |III.1| and the property stated in Proposition |III.1| We define it as follows 

T(V) := e-log [1 + (det V)™], (35) 

where k is a real positive number. Here, adj(V) denotes the adjunct of V, given by adj(V) = det(V)V“^. 

As required, the regularizing function T(V) fulhlls the following properties: 

Proposition IV.1. Let V, V' be two covariance matrices and H 6e a permutation matrix (resp., S be a symplectic 
matrix) such that V' — H^ V H fresp. V' = S'^ VS), then 


T(V') = T(V). 

Proof. See Appendix [Pj 

With the regularizing function T(V), we arrive at the following volume for sets of states: 


(36) 

□ 


Definition IV.2. Given a set of Gaussian states represented by a parameter space 0 as in Def. we define its 

volume, regularized by the functional T, to he 


Vr{V) := / T(V) vg, 
Je 


(37) 


where vq and T(V) are given by Eqs. (26),(35) respectively. 


The integral (37) is now meaningful. Indeed, as consequence of Theorem IV.I we have the following result. 


Corollary 2. The regularized volume element satisfies 

T{V)^/detg < VdetAexp(-tr[adj(V)])A™ 3 ,,,[adj(V)] 
Consequently, the integral 


log[l -I- (det vy 
(det V)™ 


(38) 


T{V)^y det gd0 


(39) 


is well-defined and bounded for any measurable subset 0 C K™ over which V is positive definite. 
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Proof. See Appendix [P} □ 

Remark IV.2. Recalling that adj(V) = det V V~^, we now have all the possible divergences in (37) suppressed: if 

^max goes to infinity then the integrand is killed to zero by the exponential ^ ); while if detV goes to zero 

log(l + x) 


then the singularity is eliminated by the logarithm log [1 + (det V)"*], using the well-known relation lim 




= 1 . 


Furthermore, thanks to Propositions III.l IV. 1 it follows that the regularized volume in (37) is invariant under 
permutation transformations and, if the states are guantum, is also invariant under symplectic transformations. 


V. EXAMPLE OF BIPARTITE STATES IN TWO-MODE SYSTEM 


We now apply the method proposed in the previous Section to a two-mode physical system {N = 2). Hence, 
the elements of the statistical model S are the pdfs in ( |T^ , where V(6) is a 4 x 4 covariance matrix and ^ = 
{qi,Pi,q 2 ,P 2 f" G r = This implies that S = {-P(C;^*)} can be at most a 10-dimensional statistical model. 

According to Remark III.1[ we consider that position and momentum variables of the same mode are not correlated. 
Thus the highest possible dimension of S reduces from to = 10 to m = 8. Then the parameter space 0 is a subset of 
the linear space M®. The Gaussian classical and quantum states are represented by 


0ci.s.ic = {9 & K®|V(0) > 0} 

= {dGK®|H(0) + zH>O}, 


(40) 

(41) 


where H is the canonical symplectic real 4x4 matrix defined in (|^. 

The separable states have to resp ect n ot just the uncertainty relation V{0) -|- > 0, but also the restriction 

V{6) -I- > 0, as stated in Theorem II.3 or equivalently V{9) -I- > 0, where H = AbH^b and Ab is the partial 


transposition defined in (14). Hence, the parameter space 0 representing Gaussian separable states is given by 


0 


separable 


= {9g K®|V(6») -f > 0, V{9) + in> 0}. 


As a consequence of Eqs. (41), (42) we also have that the Gaussian entangled states are represented by 

0 . 


^entangled ^^quantu 


- 0 , 


(42) 


(43) 


Finally, to each set of states, classical, quantum, separable and entangled, we associate a Riemannian manifold 


given by A4 = (0,G(0)), where 0 is specified by Eqs. (40), (41), (42), ( [43| ), respectively. Furthermore, G{9) is the 
Fisher-Rao metric whose components are given by (22), which is an explicit expression in terms of the covariance 
matrix V. 

However, the most general parametrization of a two-mode covariance matrix V {9) is realized through its canonical 
form and it only employs four parameters pS] , 


V{9) 


/ a 0 c 0 \ 
I 0 a 0 d I 
c 0 6 0 

\0 d 0 b J 


(44) 


where, according to our notation, the only non zero parameters are 9^ 
and = d G K. In this case, the domains of integration given by Eqs. 
the following form: 


= 0’^ = a G M, d® = = 6 G K, = c G M 

(40), (41), (42), apart from null sets, assume 


0cia,sBic = {(a, 6, c, d) G a > 0, 6 > 0, —yfab < c < -x/ofo, —yfab < d < 

0qua„tu„ = {(a, 6, c, d) G a > 1, 1 < 6 < a, -^/cf <c< y/c{, di < d < d 2 } 

U {(a, 6, c, d) G a > 1, 1 < a < 6, —y/cf < c < y/cf, di < d < d 2 } 

0=eparabia = {(o, 5, c, d) G | o > 1 , b> 1, - y/c^. < c < 0, di < d < -di} 

U {(a, 6, c, d) G K‘^1 a > 1, 6 > 1, 0 < c < yfcf,, —d 2 < d < d 2 }. 


(45) 

(46) 


(47) 
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with 


a ,, 2 ^ b , 2 1 \ I — 

-l),c,:=-(« -l).C 3 :=- - -, 

di := ■— < A :— - {oi - c^) [a6c^ - (a^ - l)((>^ - I)] , 


ab — ’ 


ab — (?■ 


Then, from the relation (221 we can compute the Fisher-Rao metric g(&) = g^vd9^ ® d9'' {g.,v = I,..., 4), whose 
components result 


5ii = 


b^{2a%^ + c‘^ + d‘^- 2ab{^ + d^)) 

2{ab — {ab — 

be bd 

^~{ab- c^Y ’ ^~{ab- d?Y ’ 


5l2 = 


922 — 


{o?Y + c^d?){c^ + d?) - Aabc^d'^ 
2 {ab — c'^Yi^b — d?Y 

Y + 

2 {ab — c'^Yio-b — d?Y 


a'^{2a?Y + d'^ — 2ab{Y + d^)) 


523 = 


ac 


{ab — YY 


, 524 = 


ad 


{ab — d?Y 


I 533 = 


ab + Y 

{ab-c^Y' 


534 = 0, 544 = 


ab + d^ 

{ab — d?Y 


At this point we are able to compute the volumes of Gaussian states whether they are classic, quantum or entangled. 
First, it is evident from (45),(461 and ( [47| ) that 0 

separable C C 0 

classic■ 

Consider now the measure V,p of Eq. (30); then, from Corollary it follows that it is absolutely continuous with 
respect to the Lebesgue measure on Also, it is monotonic. For all these reasons we have 


<d{V) UG < 


’ ©quantu 


<d{V) UG < 


HV) l^G, 


(48) 


f ©classic 


for every E e M+. Here, <?(F) = id(E — tr(F)) log [l + (det F)^] = id(E — 2(o + 6)) log I + ((a6 — c^)(a6 — d^)) 
Given the sets of states as in (|4^,(|46| and (|4^, after having obtained the set of entangled states by Eq.(|T 


we 


can also compute the volume of the latter set by the measure V 4 , of Eq.(30). In Eig. are reported the ratios of 
quantum over classical volumes, separable over classical, and entangled over classical. Figure [^clearly shows the chain 
of inclusion entangled C V,?.,separable C V<?.,quantum C V4.,classical holding true for any value of E G K+. In addition all 
the volumes go to infinity when E —>■ oo, however this takes place at higher rate for V 4 >, classical) so that all the ratios 
approach zero when E —+oo. Analogously, for E —0 all sets becomes empty and the ratios become zero. 


Ratio 



FIG. I: Ratios of volumes Vi. of Eq.(30l vs 
entangled over classical). 


E (from top to bottom quantum over classical, separable over classical, and 


Remark V. 1. The form of the covariance matrix V in (AA) is a very special one. Indeed, any 4x4 covariance matrix 
can he brought to that by a suitable transformation corresponding to some element of Sp{2,M.) x Sp{2,M.). 
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Unfortunately, the function (p just employed to get the inclusions (48) is not invariant under Sp(2,IR) x Sp(2,IR) 
transformations. 

On the other hand, the parameter space given in Eq. (42), which describes the Gaussian separable states of a 
two-mode system, has a Sp(2,K) ® Sp(2,IR) C Sp(4, K) invariant form. Hence, the function <P seems to be unsuitable 


to describe the volume of Gaussian states. In contrast, the regularizing function T of Eq.(351 is invariant under local 
symplectic transformations, as it immediately follows from Prop. IV. 1 For this reason, we propose the volume (37) 


as a suitable measure assessing differences among the Gaussian states. Again, we obtain strict inclusions among the 
volumes of classical, quantum, entangled and separable states, namely 


[ T(U) < / T(I/) ^G < / T(I/) Z.G, 

©separable ©quantum ©classic 


(49) 


and for all k G ffi+. It is worth stressing 


with T(U) = log 1 + _ c2)(a6 - (f)) 

that due to the regularization, the important quantities are not the absolute values of volumes, but rather their ratios. 

In Fig. [^are reported the ratios of the volumes of different sets of Gaussian states computed through the measure 
(37). In particular the ratios of Quantum over Glassic, Entangled over Glassic and Separable over Classic are shown 


vs At. 



FIG. 2: Ratios of volumes Vt of Eq.( 371 vs k (from top to bottom quantum over classical, entangled over classical, and separable 
over classical. 


Figure clearly shows the chain of inclusion Vt, separable C Vx, entangled C Vt .quantum C Vt, classical holding true 
for any value of k. For k —^ 0 all sets becomes empty and the ratios become zero, while for k —>■ oo they tend to 
asymptotic values. 

Both Figs, [^andj^put forward a non-monotonic behavior of the volume ratios. This effect turns out to be of purely 
geometric nature and has to be ascribed to the curved metric (by contrast one can check that with standard 
Euclidean metric the behavior is monotonic). By comparing Fig. and Fig. it is worth noticing the different 
hierarchies of volumes appearing there. This fact can be explained as follows. First we recall that a symplectic 
transformation S acts by congruence on a covariance matrix V —A SVS'^, hence the set of Gaussian states can be 
thought as the orbit for the action of Sp(2V, K) on a seed V. Then, the regularization (35) provides a “homogeneous” 


cut-off on the space of Gaussian states and as consequence a hierarchy of volumes similar to that of finite (low) 
dimensional case is obtained (in Ref. m it was shown that the volume of separable states is contained in the volume 
of entangles states, which in turn is contained in the volume of quantum states). In contrast, the regularization (29) 
provides a non “homogeneous” cut-off on the space of Gaussian states. In fact it only cuts the states arising from the 
action of elements of the non compact subgroup of Sp(2Af, K), leaving unaffected states arising from the action of the 
compact subgroup of Sp(2A^, M) (actually Sp(2A^, M) n SO(2A^, K)). Gonsequently the hierarchy of volumes obtained 
differs from that of the finite (low) dimensional case. 


VI. CONCLUSION 

In the present work, we tackled the problem of evaluating the volume of Gaussian physical states, both classical and 
quantum. The relevance of considering Gaussian states is twofold: first, Gaussian states are the most commonly 
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experimentally used CV states. Second, Gaussian quantum states are represented in the phase space picture of 
quantum mechanics as proper pdfs. Hence, Gaussian classical states are pdfs in phase space and Gaussian quantum 
states are pdfs coming from Wigner functions in phase space. Thereby, dealing with pdfs. Information Geometry 
appears as a natural and unifying approach for evaluating volume of classical and quantum states. 

By exploiting methods of Information Geometry, we associated manifolds to different sets of states; thus there 
is one manifold corresponding to classical states, one manifold to quantum states, another one to separable and 
another to entangled states. The key point in setting up such connections was that each set of states can be obtained 
by considering the pdfs parametrized by the entries of the covariance matrix. Then the manifolds are exactly the 
parameter spaces obtained by imposing constraints on the covariance matrix in order to describe classical or quantum 
states. Concerning separable and entangled states, the question is more delicate. Indeed, there is no general criterion 
to characterize multipartite entangled states. Nonetheless, by reducing to bipartite systems one could use the condition 
(191 (which turns into an operational condition for two-mode systems thanks to (HI)- Then, in this case we can also 
associate a manifold to separable states and a manifold to entangled states. Next we endowed each of these manifolds 
with a Riemannian metric, the Fisher-Rao one. Thus it was natural to define the volume of a set of Gaussian states 
as the volume (21) of the Riemannian manifold associated to it. 


Thanks to Eq. (22), we were able to show that the volume measure introduced in (21) is invariant under labelling 
permutations of modes; moreover, we proved that it has a symplectic invariant form. These results showed that the 
volume measure in (21) is suitable for estimating the volume of Gaussian states. However, since we analysed infinite 
dimensional systems, problems arose from the non-compactness of the support of the states. We overcome this 
difficulty hrst of all by resorting to an energy constraint. Hence we defined the regularizing function (29). However, 
such kind of regularization turns out to be not invariant under symplectic congruence. Then we introduced a different 
regularizing function, namely the one in (35), which came about by figuring out the functional relation given by Eq. 
. We proved that also such a function has permutation and symplectic invariant form (Prop. IV.l). 

Accordingly with these regularizing functions we have explicitly evaluated the volume of two-mode Gaussian states. 
Note that it is not the values of the volumes per se that are really relevant, but rather the ratios between the volumes 
of the various sets. As such we presented in Figs, [^andj^the ratios with respect to the volume of classical states. The 
Figures show different hierarchies of volumes with the one in Fig. [^resembling that of finite dimensional systems (at 
least for quantum states, see e.g. m)- In both cases the ratios depend on the cut-off parameter in a non-monotonic 
way due to geometric effects. This makes evident a rich structure for sets of Gaussian states. 

Finally, the presented volume measure could also be applied to three-mode systems, for which an operational 
criterion to distinguish separable states among all the quantum states is well-known |33j . Indeed, in such way the 
parameter space in Def. can be implemented for separable states as well as for classical and quantum states. 

Thereby, the volumes can be computed. Beyond that, a necessary and sufficient criterion to describe Gaussian 
separable states still lacks, hence the introduced volume measure can only be useful to provide bounds on the volume 
of sets of multipartite states. 
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Appendix A: Proof of Theorem |III.1| 


Let us notice that because of the form of P{^;9) in (15), the expression in Eq. (16) involves a Gaussian integral. 
However, before evaluating it, let us study the function 

UAO ~ d^\ogP{^-0)d,\ogP{^;9). (Al) 

By means of logarithm’s properties we can write 

2N 


log[P(e;0)]=-i 


log[(27r)^detU(d)]-b ^ 


(A2) 
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where V^if{9) is the entry a/3 of the inverse of the covariance matrix V{9). Then the derivative of Eq. (A2) reads 


5^1og[P(e;d)] = -- 


g^(det V) ^ ^ rtz-iv po 


Recall that the following relation holds 

a^(det V[9)) = det V{9) tr [V-\9) d^,{V{9))\ . 


(A3) 


(A4) 


Hence, using (A3) and (A4), we arrive at 

1 




tr[E-i(0)5^(E(0))] + ^ d^{y-p\9))Uip 

a,/3—1 
2N 

tr[V-\9)d,{V{9))]+ Y. 


(AS) 


For an analytic function /(^) and a symmetric definite-positive 2N x 2N matrix A it results 

2N 


dif{i)e 


(27r) 


2N 


det A 


exp 


1 ^ ,^-1 d d 


/l«=0: 


(A6) 


where A~-^ is the entry ij of the inverse of the matrix A and the exponential means the power series over its argument 

_^ q^ct f(Q\ 

(the differential operator). Indeed, by expanding f(^) we have that f(^) = where H = (X>/)(0), 

a\ 

cx 

with a a multi-index and T> the differential operator with respect local coordinates .. .^^v- At this point the left 
hand side of Eq. (A6) can be written as 




+S{ 


Then by performing an orthogonal transformation O, we have 

+ -\y^Dy + Y.Oy, 

where D is diagonal matrix with elements eigenvalues of A. Finally, through some algebras, we arrive at 


d^e - 


+sg_ / (27r)2^ 
det A 


which gives exactly Eq. (A6). 


Inserting the expression of P{^;9) in (15) into the relation (16) and employing Eq.(A6) we find 

1 


y(27r)2WdetE 


J =exp 


d d 




1 


Ifiu 1^—0- 


(A7) 


We are now going to evaluate the Gaussian integrals in Eq.(16) by means of the following Lemma. 

1 d d 

Lemma A.l. Let be D := - the right-hand side of (A7) we have 


9fiv{d) — ffj.,y{0) -D'^f 


(AS) 
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with 


and 


DU.k=o = -iv[v-\e)d^{v{e))\tr[v{e)d,{y-\e))\ 
+^tr [v{e)d^{v-\e))\ tr [v-\e)d,(y{e))\ 


\D^f^uk=o = \tv[v{e)d^{v-\e))]tv[v{e)d.{v-\e))] 
+^tr [v{e)d^{v-\9))v{e)d,{v-\9))\. 


Proof. From Eq. (A51, with i,j G {1,..., 2A^}, by a straightforward calculation we have 

2N 


d dL 


V d^j 


(0 = ltr[V-^d.V]+ ^ 


06,^ — ! 
2N 


,/3^i / y 

^ 2N \ / 2N ^ 

,/3=l / V/5=l / 

/ 2N \ 

+ tr [E-1 d,V] + ^ 5, (y-/) . 

\ a./3=l / 

Taking the sum over i,j G {1,..., 2A^} and evaluating the above expression at ^ = 0, we obtain 


2N 


= jYl [y ^d^v] 


, 2N 


+7 E [y-'a^y]. 


ij'=l 


Now, recall that 


2N 

1. tr [AB] = AijBij, for any pair of x iV matrices A,B; 
*7=1 


2. d^{V{9)) = 


ld9„ 


for any matrix V. 


Hence, we get Eq. (A9|. 

Eurthermore, letting i,j,h,kG {1,..., 2Ai}, we have 


d^h (aa a^z d^j ^ 

duV-^^+ a.y^fci 

Ad^V-y} d^Vr:\ 


(A9) 


(AlO) 








15 


Taking the sum over {1,..., 27V} we obtain 

E y^3Vhk + Y. d,VY d,V-^ 

\ ij,h,k ij,h,k 

+ E d.v-^ + Y d,v-^ d^vY 

i,j,h,k i,j,h,k 

+ Y + E d^V^l 

i,j^h,k i,j,h^k ) 


= i 2tr [Vd^VY tr [Vd^VY + 4tr [^ V d^VY \ ■ 


Finally, thanks to the above expression of have that the expansion in the right-hand side of Eq.(A7) 

only contains terms up to the second order. 


□ 


At this point, collecting the results in Lemma A.l together Eq. (A5) evaluated in ^ = 0, we obtain 


9^.u = \ [tr(E-i d^V)+tr{Vdf,VY] [tr (F"^ a,F) + tr (Fa^F-^)] 

+ Y (F df, F-i F d,VY ■ 


Then, the statement of Theorem IILl easily follows from relation 9^1/ ^ = —V ^{dfiV)V 


A-i 


(All) 

□ 


Appendix B: Proof of Proposition [Tll.l| 


Let us consider the permutation a : {1,...,A} —>■ {1,...,A} and the corresponding permutation matrix LI which 
entails a labelling permutation in the phase space L, i.e 11 : ^ G F —>■ ^cr G E. From the P{C,0) in (15), it follows that 
such a permutation acts on the covariance matrix F in the following manner: 


F ^ n F 


So, let V'Y s-nrl V{d) be two parametrized covariance matrices and LI a permutation matrix such that V'(0) = 
n V{6) n”. Let 0 and 0^ be the parameter spaces corresponding to V{6) and V'{d), respectively. Then there exists 
a diffeomorphism ip : Q ^ Q' with Jacobian J^p such that | det J^\ = 1. Therefore we have 

V(F') = [ de A/det5'(6») = [ dO ^/det g{9) = V(F), 

Je' Je 

where we used the equality det g'{0) = det g{9) intending g'{9) as the Fisher-Rao information matrix corresponding 
to V'{9). 

Actually, showing that detg'{9) = det g{9) we are proving a stronger relation between the Fisher-Rao metrics 
corresponding to F and F^ In fact the following relation holds true, 

tr [V dp {V'Y y' d^{V')Y = F [nFR"^ dp (nFn'^)-^ RFH^ d^{UVXF)Y 

= tr [F dp F-i F d^VY . 


where we used the independence of 11 from 9 and the invariance of the trace under cyclic permutation. 

Then, recalling the relation (22), we arrive at = gp,j, for every g,,v G {l,...,m}. Here, denotes the 
component of the metric corresponding to V'. Thereby, detg'(0) = detg(0) trivially holds true. 


Focusing on the quantum states, it is well-known that the uncertainty relation F -I- > 0 is invariant under 

symplectic transformation [25]. So, let us consider two parametrized covariance matrices V'(9) and F(0) and a 
symplectic matrix S such that ViS) = S V{9) S'^. Then, the parameter spaces and 0q„antum: corresponding 

to those different matrices, coincides. 
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Furthermore, we have 

tr [F' {V')-^ V d^iV')-^] = tr {SVS'^)-^ SVS^ d^{SVS^)-^] 

= tr [y v-^ y d,v-^], 


where we used the independence of S from 9 and the invariance of the trace under cyclic permutation. Thus, from 
(22) it trivially follows that, 

V(y') = f dO ^detg'iO) = f d9 i/det g{9) = V(y), 

©_..._ 0riiianeii7-n 


where g {6) and g{6) denote the Fisher-Rao information matrix corresponding to the covariance matrices V' and V, 
respectively. □ 


Appendix C: Proof of Theorem |IV.l| and Corollary [T] 


The proof of Theorem IV. 1 relies on an inequality of the determinants of Gram matrices [34j . We recall that for an 
inner product space ?{ and a set of points {cci,... ,Xn} in H, the Gram matrix of this set is the n x n matrix defined 
by G{xi,...,Xn) = 


Lemma C.l. Let "Hi,7^2 be two inner product spaces and T : "Hi —> 772 be a bounded linear operator. Let 
{xi,... ,Xn} € 77i be an arbitrary set. Then 


AetG{Txi,. .. ,Txn) < ||T|p” det G(xi,..., a;„), 


(Cl) 


where ||T|| denote the operator norm of T. 

Proof of Theorem \IV. 1\ Since V is symmetric, positive definite, its adjunct, given by adj(y) = det(y)y“^ is also 
symmetric, positive definite. We have 

ffu!' = = ^ ^^^^^^^^ tr[adj(y)((9^y)adj(y)(a^y)] 

1 

“ (dety)2^'""’ 
where the matrix is given by 

ffw' = ^tr[adj(y)(a^y)adj(y)((9,,y)] 

= itr[(adj(y)C2(a^V)adj(y)C2)(adj(y)i/2(a^y)adj(y)i/2)] 

= i((adj(y)C2(a^V)adj(y)C2), (adj(vy/^(d.v)adj(vy/^)}^, 

with (•,-)f denoting the Frobenius inner product given by {A,B)p = tr [(i7^A)] for any matrices A,B with same 
dimension. 

Consider the linear operator Ty : M 2 Afx 2 Af ^ 2 Ny. 2 N^ with R 2 JVx 2 Af ^j^g Frobenius inner product, defined 

by 

TyA = adj(y)C2Aadj(y)C2. (C2) 


Then 


IITvAIIf = ||adj(y)C2Aadj(y)C2||^ < ||adj(y)C2|| ||Aadj(y)C2||^ 

< ||adj(y)|| ||A||f, 


with equality if A = 7, where we have used the property that ||adj(y)|| = ||adj(y)^/^|p by the symmetric, positive 
definiteness of adj(y). Thus 


||Ty|| = ||adj(y)||. 


(C3) 
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Then we have 


~g^. = -{T{d^Vind,V))p. 


Let E be the m x m matrix defined by 


By Lemma [C.l[ 


It follows that 


= \tT[{d^V){d,V)] = \{{d^V),{d,V))F. 


detg< ||Ty|p™detL;= ||adj(y)|p™det£;. 


, ||adj(F)||2™ , ^ 

det g < E ————— det E = 

(dety)2”^ 


niadj(y)|| \ 
I det 1/ / 


2m 


det L; = f det E. 

\ det 1/ ) 


Let {Xk}k=i be the eigenvalues of V. From the relation adj(y) = det(F)t7 it follows that the eigenvalues of adj(y) 

f det(y) ^ , 

< —--> and thus 

I J k=i 


are 


Amax(adj(y)) = ^ det5 < 


Xn^iniV) 


2m 


det E. 


□ 


This completes the proof of the theorem. 

Proof of Corollary^ The first expression of the Corollary follows from the bound given in Theorem IV. 1 and the 
definition of 

We now show that the integral 

I 'P{V)\/det gd9, 


log(l -|- x\ 

is bounded. By the inequality log(l + a:) < a; for all a: > 0 and the limit lim- = 1, we always have 

x—^0 X 


log[l + (det VY' 
(det V)"* 


< 1 whenever det V > 0. 


Consider now the factor iJ(E — tr(V))A™^x(adj(C)). As in the proof of Theorem IV.1 let {\k}k=i be the eigenvalues 
of V, arranged in decreasing order. Then 


Amax(adj(V)) = 


det V 
X 2 N 


2N-1 


n 

i=i 


'^2Af-r > \ 2JV-r 

2V- 1 / 


< 


tr(V) 
2V- 1 


27V-r 


< 


/ E 


V2A^- 1 


2A7-r 


where E is a suitable real positive constant bounding from above tr(V), which is a positive linear function in 9. 
Furthermore, the domain of integration is now bounded from above, indeed it is given by {0 G | V(0) > 0} C 
{0 G K™ I tr(V) < E}. As sort of evidence, let us consider the following expression for det(V) [55] , 


,_, 7_iVi+i 

det(A^)= E 


fci. .,k/^ 1—1 


(C4) 


N 


where the sum is taken over the set of all integers fc; > 0 satisfying the equation iki = N. Now, whenever V is 


i=i 

positive definite matrix, through Cauchy-Schwarz inequality we have that t r(V^ ) < tr(V)^; thus, by induction on n, 
with 1 < n < A^, we arrive at tr(V*) < tr(V)* for every I < N. So, from Eq.(|C4[) we have that. 


N 


0<det(V)< ^ 


(C5) 

/.'"t k:i ! 

ki ,,fciv ^ = 1 

whenever V is positive definite. Moreover, because of trace of a principal minor of V is smaller than tr(V) then, we 
earn same bounds as in (C5) for determinant of every principal minors whenever V is positive definite. Thus, we have 
the convergence of the integral. □ 
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Appendix D: Proof of Proposition |IV.1| and Corollary 

Proof of Proposition IV. 1 Consider a permutation matrix 11 such that V' = I\^W. Then, because of the 
property of unitary determinant det If = det 11^ = 1, and the fact that the permutation matrices If do not depend on 
the parameters we have 

det V' = det (n^FIl) = det 11^ det V det 11 
= det V, 

tr [(detP')(^')”^] = detP' tr [(P')”^] 

= tr [(detP)(P)-i] . 


Hence, from Eq. (35l it immediately follows that T(E(d)) = T(E'(d)). 

In the same way, since det ^ = 1 = det S^, with S a symplectic matrix such that V' = S^VS, then we have 

det V' = det V, 

tr [(detE')(^')”^] = tr [(detE)(H)-i] . 

Hence, from Eq. ([^ it immediately follows that T(E(0)) = T(E'(0)). 


□ 

Proof of Corollary^ The first expression of the Corollary follows from the bound given in Theorem |IV. 1| and the 
definition of T. 

We now show that the integral 


T {V)^/det gd0 


log(l -|- x) 

is bounded. By the inequality log(l + a:) < a; for all a: > 0 and the limit lim- = 1, we always have 

x-fO X 


log[l + (det VY' 
(detF)"* 


< 1 whenever det E > 0. 


Consider now the factor exp(—tr[adj(E)j)A™ 3 ^,j(adj(E)). As in the proof of Theorem IV.1 let be the eigen¬ 

values of V, arranged in decreasing order. From the relation adj(E) = det(E)E“^, it follows that the eigenvalues of 


2N 


adj(E) are | \ ^ have 

I J k=i 


tr[adj(E)] = ^ 


2N T , 2iV 2N 

det V 






“ X! n 

k=l j = l,j^k 


which is a positive polynomial in the parameters Furthermore 


Amax(^dj(E)) — 


detV 

^2N 


2N-1 


E 2N —1 \ ' 

2iV- I 


n - 


< 


i=i 


2iV-l 


< 


tr(V) 

2iV- 1 


2N-1 


where tr(y) is a positive linear function in 6 . Thus as 0 grows, the expression 

exp(-tr[adj(E)])A™,,^(adj(E)) 

decays exponentially, leading to the convergence of the integral. 


□ 


111 K.S. Mallesh, S. Chaturvedi, R. Simon, N. Mukunda, States of Physical Systems in Classical and Quantum Mechanics. 
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